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# I /usr/local /bin/perl 

# COPYRIGHT (c) 1598 Telcordia Technologies Inc., 

# All Rights Reserved. 
# 

# PROPRIETARY - BELLCORE AND AUTHORIZED CLIENTS ONLY. 

# This document contains proprietary information that shall 

# be distributed or routed only within Telcordia Technologies 

# (Telcordia)/ and its authorized clients/ except 

# with written permission of Telcordia. 
# 

n $Id: getHostLoc.pl,v l.X 1999/05/20 22:27:07 rmartija Exp rmartija $ 



require 'getopts.pl' ; 
undef ; 

SUSAGE = "Usage: " . $0 . " I-D] -u file -m file\n" . 
"Flags ;\n" , 

" -D debug mode\n" . 

" -u file file containing the list of unclassified IP\n" . 
It addresses (i.e, those with un3cnown locations) \n" , 

» and their characteristics. \n" . 

" -fti file file containing the means and inverse of covariance\n" 
" matricesXn" . 

"Examples :\n^ . _ . „ _ . 

" SO -u unknowns -m matrix" ; 



%g_means = ( ) ; 
%g_inverse = ( } 
®g_locales; 
$g_debug; 
$g_at tributes; 



# 

# 

siib getDistance { 

TTiy( Sloe, Sdata ) = 

my( ®X ) = ®$data; 

my( ®mu ) - ©{$g_roeana{$loc} } ; 

my( ©Sigma ) = {e){$g_inverse{$loc}} ; 

my ( Qdiff , ©prod ) ; 
my ( $i/ $ j ) ; 

for( $i ^ 0; $i <" Sg_attributes; $i++ ) { 
$diff[$i] = $mu[$i] - $Xt$i]; 

} 
# 

# compute diff (transpose) * sigma. diff (transpose) is a 1 x N matrix 

# and Sigma is a N x N matrix, the result is a 1 x H matrix. 



PAGE 28«3* RCVD AT imm 10:50:22 AM [Eastern Daylight T^^ 



OCT-23-04 10= 14 FROM : TELCORDIA LEGAL DEPT. 



ID: 17323303004 



PAGE 23/33 



# 

tori $i = 0; $i <a $g_attribut:es; $i++ ) { 
$prod[$i] = 0. ; 

for( $j = 0; $j <^ $g_accribures; $j++ ) { 
$prod[$i] += $diff[$j] * $sxgTna [$i] / 

} 

) 

# 

^ multiply the matrix obtained above, l.e prod, wich diff. prod is a 
# a 1 X N matrix and diff is a N x 1 matrix, the result is a scalar, 
# 

ticy( $dist ) =0; 

foci $1 = 0; $i <- $g__attributes; $i++ ) { 
$dist 5prod[$ir * $diff[$i]; 

) 

recurn $di3t; 

} 



# - 

# 

9uh readMeansAndMatrices { 
my( $file ) ^ 

open(-F. "< $file" ) ; . ^ _ _ . 

©lines - <F>; 
olose ( P ) ; 

my( $n_row6, $cur_row, $line^num ) = (-1^ 0, 0); 
my( $cur_loc, $n_mean6 ); 

foreach( ©lines ) { 
chop; 

$lliie_num++; 

next if /^\s*$/; # skip blank lines 

if ( $_ =~ /"US.*:\s.*(-*)/ ) { 

die "ERROR: $file is corrupted\n-> line $line_numj $_\n" 
unless $n._rows < 0; 

# $1 contains the state string (e.g, NJ) 
$cur_loc = "$1,US"; 

$ciir_row = 0; 

} 

elsif ( $_ /*NOHUS.*:\s*(.*)/ ) { 

die "ERROR: $file is corrupted\n- > line $line_num: $_\n" 
unless $n_^rowB < 0; 

# $1 contains the country string (e.g. BE) 
$cur_Xoc ^ •'$1,$1"; 

$cur_row = 0; 

elsif ( $_ /"MEAN.*:\s*(.*)/ ) { 
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die "ERROR: $file is corrupted\a-> line $line_auro: $_\n" 
unless $n_row9 < 0; 

# $1 contains something like 18.43 1X30,71 20-00 170-71 19.57 228.5 
Thy( ©means ) = split { ' $1 ) ; 

$n_means = $n^rov;3 = $#means; 
$g_means{$cur_loc} « \©mean3; 

elsif ( =^ /-^INVERSE. »:\6*(.*)/ ) { 

die "ERROR: $file is corrupted\n-> line $line_num: $_\n" 
unless $cur_row == 0; 

elsif ( S_ =" /"([A-za-z]+) .*:/ ) { 

die "ERROR: Invalid Tag in $file\n-> line $line_num: $_\n"; 

else { 

Tny( ©row ) = aplic( •'.$_)/ 

# make sure the matrix Is a $n_means X $n_meana array 
die "ERROR: $file is corrupted\n-> line $line_num: $_\n" 

unless $#row $n_mean3 $cur_row <:= $n_means; 

Tny( $r_Bntry ) = (®row] ; 

push( ®{$g_inverse{$cur_loc}} , $r_entry ); 



} 



$cur_rowi"*-; 
$n rows--; 



} 

die "ERROR: $file is corrupted. More data expected. \n" unless $n_row3 < 0; 

©g^looalea - keys %g_mean3f 
return $n_means/ 



# - 

# 

sub classifylPs { 
my( Sfile > ^ 

open( F, $file» ) ; 

my( ®data, $tloC; $loc, %dist, Smin ); 



while ( <F> ) { 

chop; 

next unless S_ =- (\d+) \ . (\d*) \ . {\d-h) \ . (\d+) , *\s* ( . *) /; 
($ip, ©data) - split( ' ' ); 
next unless §#data == $g_attribute3/ 

Smin = time; ft initialize $dist to some arbitrary large number 
# such as the number of seconds since 1/1/1970 
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foreach $tloc ( ®g_locales ) { 

Sdist{$tloc} ="&getDistance ( $tloc, \@data ); 
i£( .$dist { $tloc} < $min ) { 
$min - $di6t{$tloc}/ 

$loC = $tlOC; 



if ( $g.debug ) { 

foreach $key (sort keys %dist) { 

printt "%-15s %-8s %7.2f\n», $ip, $key, $di3t{$key}; 

} ' 

} 

printf "%-15s %-8s\n", $ip, $loc; 

} 

Close ( P ) ; 

} 



####################### main program #########»################## 
################################################################# 

$x = &Oetopts( 'uimiD' )/ 

die "$usAGE\ii" -unless ($x ne 

die "$aSAGE\a" unless ($opt_u && $opt_tn) ; 

die "ERROR: cannot open $opt_u\n" unless -e $opt_u; 
die "ERROR; cannot open $opt_m\n" unless -e $opt_m; 

$g_^debag = 1 if ( $opt_D ) ; 

$gj[^attributes = sreadMeansAndMat rices ( $opt_m ) ; 
tclassifylPs ( $opt_u ) ; 
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